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The  usual  procedure  for  calculating  inhomogeneously  broadened  line 
shapes  of  rare-earth  ions  in  crystalline  solids  is  to  assume  a  par¬ 
ticular  form  of  the  interaction  Hamiltonian  and  to  proceed  by  using 
perturbation  theory. *  The  procedure  outlined  here,  although  a  form  of 
perturbation  theory,  is  presented  in  a  form  that  will  allow  for  a 
simpler  method  of  ascertaining  the  magnitude  of  various  forms  of  inter¬ 
actions.  This  method  is  particularly  applicable  if  computational  facil¬ 
ities  are  available  to  calculate  the  crystal-field  energy  levels  of 
rare-earth  ions  in  various  host  crystals. 

The  analysis  of  rare-earth  ions  in  solids  assumes2  that  the  Hamil¬ 
tonian  is  composed  of  two  parts,  HQ  and  H^ •  The  first  part,  H^ ,  con¬ 
tains  the  "free  ion"  parameters;  these  are  chosen  so  that  in  the  absence 
of  the  crystal  field  an  adequate  representation  of  the  centroids  of  the 
various  multiplets  is  obtained.  The  second  part,  ,  is  the  Hamiltonian 
for  the  crystal  field,  written 


H1 


^  BkqCkq 


(1) 


where  the  crystal-field  parameters,  ,  are  such  that  k  is  even  and 
-k  <  q  <  k.  A  sun  over  the  electrons  of  the  configuration  4fN  is 
implied  in  equation  (1)  but  will  be  suppressed,  except  when  necessary 
for  clarity.  The  represent  the  effect  of  the  surrounding  host  on 

the  rare-earth  ion.  It  is  generally  assumed  that  Hq  is  not  sensitive  to 
local  variations  within  a  crystal  and  that  the  individual  Stark-split 
rare-earth  energy  levels  are  most  sensitive  to  changes  in  the  B^g.  He 
shall  assurve  so  here.  = 


For  a  small  change  in  the  ,  we  can  then  write  for  the  shifted 
energy  level 


rl 


.^n 


1B.  DiBartolo,  Optical  Interactions  in  9olids,  John  Wile y  and  Sons , 
Inc.,  New  York  (1968), 


2B.  G.  Wybourne,  Spectroscopic  Properties  of  Rare  Earths,  Wiley -on/ 
Interscience,  New  York  (1965). 
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where  H|n>  =  E  |n>,  H  *  Hn  +  Hi  ,  and  the  superscript  indicates  the 
n 

evaluation  of  the  quantities  in  the  absence  of  any  change  in  B^q.  The 
change  in  the  energy  separation  between  two  energy  levels  n  and  m  due  to 
a  change  in  the  can  be  written 
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Now  if  we  perform  an  average  appropriate  for  a  particular  physical 

model  of  the  6b  ,  we  obtain  for  the  line  shift 
kq 
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nm 


nm 


0  v-' 

-  w  =  V 
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kq 


C  6Bv 
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(4) 


where  the  f”  are  independent  of  the  particular  ensemble  of 
Similarly,  averaging  the  square  of  equation  (3)  results  in 


nm.nm 
kq  k’q* 


6b  6b 

kq°  k’q' 


(5) 


Squaring  the  result  in  (4),  subtracting  from  (5),  and  reducing  the 
resulting  expression,  we  have 
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*  v  '  nm  nm  / - - \ 

r2  =  )  f  f  ( 6b  6b  -  6b  6b  ) 
nm  L->  kq  k'q'  \  kq  k'q*  kq  k'q '  / 


where  T2  =  u>2  -  ( )  .  The  quantity  T2  is  the  second  moment  of  the 

nm  nm  \  nm/  1  nm 

distribution  of  shifts  and  can  be  a  useful  measure  of  the  line  width. 
Higher  moments  of  the  distribution  may  be  calculated  via  the  above 
procedure • 


The  result  given  in  equation  (6)  for  the  line  broadening  is  very 
general  at  this  point.  The  resultant  broadening  theory  can  be  applied 
to  impurities,  concentration  effects,  lattice  vibrations,  and  a  variety 
of  other  effects.  The  particular  type  of  broadening  determines  the 
origin  of  the  6b^  and  the  type  of  averaging  to  be  considered  appropri¬ 
ate  for  that  particular  mechanism.  The  derivatives  3En/3B^  which  are 
common  to  any  mechanism  can  be  readily  calculated  by  using  Feynman' s 
theorem,3  together  with  the  wavef unction  characterizing  the  level  ,  so 
that  the  calculation  of  the  f^  is  straightforward. 

As  a  specific  example  of  the  use  of  the  results  given  in  equations 

(4)  and  (6),  we  consider  an  impurity  rare-earth  ion  (or  one  of  the 

constituents)  in  a  solid  in  which  a  rare-earth  ion  is  a  stoichiometric 

constituent.  Consider  an  electron  at  r  on  the  impurity  ion  at  the 

♦ 

origin  and  an  electron  at  r^  from  the  nucleus  of  a  host  rare-earth  ion 

at  R  .  The  interaction  of  these  two  electrons  can  be  written4 
P 

H'  *  — - - - - —  *  e2  £  (-1)*  M  <k(q)a(a)  |k+a(q+a)> 

I R  -  r  +  rj  |  kq  L  J 


X  rk<W*)r*Caa(?l) 


_k+a+ 1 


3B,  Merzbacher,  Quantum  Mechanics ,  John  Wiley  and  Sons,  Inc,,  New  York 
(1961),  p  396. 

4S.  R.  Judd,  Angular  Momentum  Theory  for  Diatomic  Molecules,  Academic 
Press,  New  York  (197 5),  p  98. 


where  (^^k8)  a  bi-nom*-al  coefficient,  <k(q)a(a) |k+a(q+a)>  is  a 

Clebsch-Gordon  coefficient  and,  as  in  equation  (1),  -a  <  a  <  a.  Also, 
as  in  equation  (1)  the  sum  over  the  electrons  of  both  rare-earth  ions  is 
assumed,  'Bie  result  given  in  equation  (7)  can  be  written 


so  that  the  total  due  to  all  the  ions  is,  summing  over  host  sites 

P# 


where  we  have  used  the  fact  that  the  Clebsch-Gordon  coefficient  in 

equation  (7)  can  be  written  in  terms  of  binomial  coefficients. 5  If  we 

further  assume,  as  in  the  point  charge  model,  that  the  B.  are  related 

k  K 

to  the  K,„,  we  can  write  6B,  -  <x  >  6A.  where  the  <r  >  contain  the 

kq  kq 

shielding  factors  and  other  quantities  used  in  the  point  charge  model. 


The  rare-earth  ions  of  the  host  material  are  in  particular  energy 
levels,  so  that  we  can  replace  r*C#o(rj)by  its  expectation  value  as 

^.fp)  ■  <)lr*c..(ii)u>  •  "01 

where  the  subscript  p  on  j  indicates  the  rare-earth  ion  at  Rp*  The 

quantity  -eQ  (j)  is  the  multipolar  moment  of  the  host  rare-earth  ion  in 
SCI  _ 

the  electronic  state  j.  Let  Q  be  the  thermal  average  of  the  multi- 

ao 

polar  moment  at  site  p. 


■  Z  •Xp(-60E)p)  °.a(5p)/z  ■ 

5m.  e,  Rose,  Elementary  Theory  of  Angular  Momentum,  John  Wiley  and 
Sons,  Inc,,  Net*  York  (1957), 
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where  Z  =  £  exp  and  0q  *  1/kT.  The  average  Q^a  is  independent 

of  p  if  we  neglect  the  interaction  of  the  constituent  ions.  Also  in 
equation  (11),  degeneracies,  if  they  exist,  must  be  taken  into  account 
in  both  Q  and  Z.  Using  in  equation  (9)  the  result  given  in  (11),  we 
can  write 


52  Gka  Q+  S. 

*-*  qo  aa  k+a,q+a 
aa 


where 


qo 


(-1)' 


(12) 


and 


NM 
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Since  dB,  *  <r  >  <5a,  ,  we  can  use  equation  (12)  in  (4)  to  obtain  the 
kq  kq 

shift  as  a  function  of  temperature  in  the  transition  of  the 

rare-earth  ion  at  the  origin.  This  shift  contains  a  term  that  persists 
at  absolute  zero  which  arises  from  the  fact  that  the  multipolar  contri¬ 
butions  have  not  been  contained  in  the  zero-order  crystal- field  parame¬ 
ters.  This  shift  at  absolute  zero  can  be  subtracted  out  so  that  the 
relative  shift  is  meaningful. 


In  obtaining  the  average  of  , »  terms  of  the  form 

®aa(^p)®b0 Op')  in  'the  cross  products  occurring  in  equa¬ 

tion  (9).  The  average  values  of  these  products  can  readily  be  shown  to 
be 


®aa(?p)®bB@p')  *  ®aa®b{$  ^pp'  +  ~  ^pp')  ®aa  ®bB  ' 
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where  the  first  terra  is  the  thermal  average  of  a  product  of  multipolar 
moments  on  the  same  site.  It  is  assumed  that  the  multipolar  moments  at 
different  sites  are  uncorrelated.  Using  the  result  (13)  we  have 


*W\ 


'q' 


-  6k,  6k, 


kq 


aa 
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(14) 


The  result  (14)  may  then  be  substituted  into  (6)  to  obtain  T _ .  The 

run 

summations  occurring  in  equations  (12)  and  (14)  are  similar  to  those 
arising  in  the  point  charge  model;  the  number  of  terms  is  restricted  by 
the  point  symmetry  of  the  site  occupied  by  the  rare-earth  ion.  The 
expectation  values  in  equation  (10)  can  be  evaluated  after  the  Hamilton¬ 
ian  for  the  host  rare-earth  ion  has  been  diagonalized.  By  the  use  of 
Feynman's  theorem,  the  Q  are  given  by 


cnOj> 


<ra> 


(15) 


where  the  Bao  are  the  crystal-field  parameters  for  the  host  rare-earth 
ion. 
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